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Abstract 

Background: Several genome-wide association studies (GWAS) involving European populations have successfully identified 
risk genetic variants associated with type 2 diabetes mellitus (T2DM). However, the effects conferred by these variants in 
Han Chinese population have not yet been fully elucidated. 

Methods:\Ne analyzed the effects of 24 risk genetic variants with reported associations from European GWAS in 3,040 Han 
Chinese subjects in Taiwan (including 1,520 T2DM cases and 1,520 controls). The discriminative power of the prediction 
models with and without genotype scores was compared. We further meta-analyzed the association of these variants with 
T2DM by pooling all candidate-gene association studies conducted in Han Chinese. 

Results: Five risk variants in IGF2BP2 (rs4402960, rs1470579), CDKAL1 (rs10946398), SLC30A8 (rs13266634), and HHEX 
(rsl 1 1 1875) genes were nominally associated with T2DM in our samples. The odds ratio was 2.22 (95% confidence interval, 
1.81-2.73, P<0.0001) for subjects with the highest genetic score quartile (score>34) as compared with subjects with the 
lowest quartile (score<29). The incoporation of genotype score into the predictive model increased the C-statistics from 
0.627 to 0.657 (P<0.0001). These estimates are very close to those observed in European populations. Gene-environment 
interaction analysis showed a significant interaction between rsl 3266634 in SLC30A8 gene and age on T2DM risk 
(P<0.0001). Further meta-analysis pooling 20 studies in Han Chinese confirmed the association of 10 genetic variants in 
ICF2BP2, CDKAL1, JAZF1, SCL30A8, HHEX, TCF7L2, EXT2, and FTO genes with T2DM. The effect sizes conferred by these risk 
variants in Han Chinese were similar to those observed in Europeans but the allele frequencies differ substantially between 
two populations. 

Conclusion: We confirmed the association of 10 variants identified by European GWAS with T2DM in Han Chinese 
population. The incorporation of genotype scores into the prediction model led to a small but significant improvement in 
T2DM prediction. 
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Introduction 

Type 2 diabetes mellitus (T2DM) is a complex disease 
influenced by both genetic and environmental factors. The 
heritability of T2DM is relatively strong with an estimated h 2 of 
31-69% [1]. Previous genetic studies have suggested the 
involvement of multiple genes with modest effects in the 
pathogenesis of T2DM [2] . This notion was supported by several 
genome-wide association studies (GWAS) for T2DM in European 
population [3-7] . These GWAS showed associations of approx- 
imately ~40 risk variants with T2DM in European population. 
Further large-scaled meta-analyses confirmed these associations 



and estimated their relative contributions in European populations 
[8] . These discoveries greatly advanced our understanding toward 
the genetic architecture of T2DM and provided valuable tools for 
the prediction of personal T2DM risk in European populations. 

The prevalence of diabetes mellitus has increased rapidly in 
Chinese populations in recent decades. In 2013, the prevalence of 
diabetes and prediabetes was estimated to be 11.6% and 50.1% 
respectively, suggesting that there were 113.9 million Chinese 
adults with diabetes and 493.4 million with prediabetes [9]. This 
dramatic surge of T2DM prevalence poses a serious threat to the 
public health of Chinese populations. Since diabetes can be 
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effectively prevented by life-style or pharmacological intervention 
in high-risk patients, it is important to identify high-risk subjects 
for preventive measures. With the strong heritability of diabetes, 
genetic information is expected to offer additional benefits towards 
the identification of high-risk subjects. Previous studies incorpo- 
rating genetic scores into T2DM prediction models have 
successfully demonstrated the benefit of utilizing such approach 
[1 0, 1 1] . However, this approach has not yet been validated in Han 
Chinese population. Given the heterogeneous genetic structures 
between European and Chinese populations, it is essential to 
confirm the association and the predictive value of these genetic 
variants in the Chinese. 

In this case-control study, we genotyped 24 risk variants 
identified from European GWAS in 3,040 Han Chinese. The 
associations of these variants with T2DM were analyzed in our 
sample and were further validated by a meta-analysis pooling 20 
case-control association studies of Han Chinese. The discrimina- 
tive power of the prediction models with and without genotype 
information were then compared. 

Materials and Methods 

Study participants 

A total of 760 T2DM patients were recruited from the 
metabolism clinics of the National Taiwan University Hospital 
(NTUH) and another 760 T2DM patients were recruited from the 
metabolism clinic of the Yunlin branch of NTUH. T2DM were 
diagnosed according to the criteria of the American Diabetes 
Association [12] or the use of anti-diabetic therapy. Patients with 
ages of onset below 35 years were excluded. In addition, 760 
glucose-tolerant healthy controls were recruited from the health 
check-up service of NTUH and another 760 controls were 
recruited from a community screening for metabolic syndrome in 
the Yunlin county of Taiwan. Glucose tolerance was defined as 
fasting plasma glucose < 126 mg/dl or 2-hour plasma glucose < 
200 mg/dl during a 75-g oral glucose tolerance test (OGTT). 
Written informed consent was obtained from each participating 
subject, and the study was approved by the institutional review 
board of the National Taiwan University Hospital. 

Selection of SNPs and genotyping 

Twenty-four genetic risk variants were selected from GWAS or 
well-established candidate-gene association studies for T2DM in 
European populations [3-7]. In view of the low risk allele 
frequencies and the negative T2DM association of rs7903416 in 
the TCF7L2 gene observed in previous researches in Chinese 
population [13,14], we genotyped another SNP in this gene, 
rs290487, which has been reported to be associated with T2DM in 
Chinese [13]. Genotype data of rs7903146 were retrieved from 
our previous study [13]. Genotyping was performed using the 
GenomeLab SNPstream genotyping platform (Beckman Coulter) 
and its accompanying SNPstream software suite. The concordance 
rate based on this platform was 99.62% [15]. 

Search criteria for meta-analysis 

We searched the PubMed database with the keywords 
"genetic", "association", "genome", "genome-wide association 
study", "type 2 diabetes mellitus", "type 2 diabetes", "diabetes", 
"Chinese", "Taiwanese", and "Han Chinese" before October 
2013. Eligible criteria included (1) candidate-gene association 
studies for T2DM published in peer-reviewed journal, (2) 
participants recruited from independent populations, (3) genotype 
information sufficient to calculate allele counts, (4) candidate genes 
tested in our study and (4) study population restricted to Han 



Chinese. A total of 21 studies that fulfill were indentified [13,16- 
35] . Full-text manuscripts were obtained for all studies. One study 
was excluded because its study population was confined to She 
Chinese population [35]. Data extraction was performed and 
examined independently by two reviewers (YCC and YHY).A11 
procedures were conformed to the Preferred Reporting Items for 
Systemic Review and Meta-Analysis (PRISMA) guidelines for 
meta-analysis [36] . The PRISMA check list and flow diagram was 
shown in Checklist SI. 

Statistical analyses 

Hardy-Weinberg equilibrium (HWE) test was performed for 
each SNP in the control group before marker-trait association 
analysis. Tests for the associations of each SNP and haplotype with 
type 2 diabetes were conducted using logistic regression. Nominal 
two-sided P-values were reported. Multivariate analysis with age, 
gender, and BMI as covariates was performed using multivariate 
logistic regression. The odds ratio (OR) and 95% confidence 
interval (CI) associated with each risk allele were also estimated. 
Pairwise gene-gene and gene-environment interactions were 
analyzed by logistic regression. The significance of interaction 
was adjusted for multiple testing using the Bonferroni method. 

To test the cumulative effects of genetic variants on the T2DM 
risk, weighted genetic score for each risk allele was calculated using 
the beta-coefficients of logistic regression model. AH participants 
were divided into four equal groups according to their genetic 
score (<29, 29-31, 31-34, >34). The OR and 95% CI for each 
group were estimated using the lower quartile group (score < 29) as 
the reference group. The statistical power of this study for each 
SNP was estimated using the Genetic Power Calculators (http:// 
pngu.mgh.harvard.edu/~purcell/gpc/) assuming diabetes preva- 
lence of 8% [37]. Meta-analysis under fixed effect models were 
used to estimate pooled odds ratio (OR) using the Comprehensive 
Meta-Analysis software (Biostat, Englewood, NJ). Cochran's 0_test 
and / was used to assess heterogeneity between the individual 
studies. The Z test was used to determine the significance of the 
pooled OR. 

Results 

Characteristics of study subjects and SNP information 

Twenty four SNPs were successfully genotyped in 1,502 
unrelated T2DM cases and 1,518 glucose-tolerant controls except 
for rsl081 1661 in the CKD2A/B gene, in which genotyping failed 
in all samples. The demographic and biochemical characteristics 
of the study participants are shown in Table 1. Basic information 
of these SNPs is summarized in Table 2. All SNPs were in Hardy- 



Table 1. Characteristics of study participants. 





Case 


Control 


P value 


Number 


1502 


1518 




Age (year) 


60.42 ±11. 83 


55.83± 15.83 


<0.0001 


Male sex (%) 


51.26 


50.86 


0.82 


Body mass index (kg/m 2 ) 


25 .45 ±4.2 7 


24.27 ±3.66 


<0.0001 


Fasting glucose (mg/dL) 


15 7.00 ±54.87 


93.84± 17.27 


< 0.0001 


Triglycerides (mg/dL) 


17 1.82 ±88.85 


115.76±77.99 


<0.0001 


Total cholesterol (mg/dL) 


195.02±50.13 


197.18±38.15 


0.20 


doi:1 0.1 371 /joumal.pone.0095045.t001 
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Table 2. Information of SNP. 





SNP rs# 


Gene 


Chr. 


Chr. position 


Gene region 


Allele 




Risk 
allele 


Controls HW 

P 


Call Rate 

(%) 














Minor 


Major 








1 


rs2641348 


ADAM30 


1 


120149926 


Exon L [Leu] => P[Pro] 


C 


T 


C 


1 


99.5 


2 


rs 10923931 


NOTCH2 


1 


120319472 


intron 


T 


G 


T 


1 


99.7 


3 


rs7578597 


THADA 


2 


43644474 


Exon T [Thr] ^ A[Ala] 


C 


T 


T 


1 


99.9 


4 


rs 10490072 


BCL11A 


2 


60581582 


intergenic 


C 


T 


T 


1 


99.1 


5 


rs4402960 


IGF2SP2 


3 


186994389 


intron 


T 


G 


T 


0.45 


98.4 


6 


rs1470579 


IGF2SP2 


3 


187011782 


intron 


C 


A 


C 


0.65 


98.6 


7 


rs17036101 


SYN2.PPARG 


3 


12252845 


intergenic 


A 


G 


G 


0.64 


99.6 


8 


rs4607103 


ADAMTS9 


3 


64686934 


upstream 


T 


C 


C 


0.91 


99.1 


9 


rs10010131 


WFS1 


4 


6410977 


intron 


A 


G 


G 


1 


99.4 


10 


rs6446482 


WFS1 


4 


6413755 


intron 


C 


G 


G 


0.54 


99.4 


11 


rs 10946398 


CDKAL1 


6 


20769003 


intron 


C 


A 


C 


0.42 


99 


12 


rs9472138 


VEGFA 


6 


43919730 


intergenic 


T 


C 


C 


0.74 


99 


13 


rs864745 


JAZF1 


7 


27953796 


intron 


G 


A 


A 


0.23 


98.6 


14 


rs 13266634 


SLC30A8 


8 


118253964 


Exon R [Arg] =M/V[Trp] 


T 


C 


C 


0.35 


99.2 


15 


rsll 1 1 875 


HHEX 


10 


94452862 


downstream 


C 


T 


C 


0.80 


99.3 


16 


rs7923837 


HHEX 


10 


94471897 


intergenic 


G 


A 


G 


0.79 


98.7 


17 


rs7903146 


TCF7L2 


10 


114748340 


intron 


T 


C 


T 






18 


rs290487 


TCF7L2 


10 


114899720 


intron 


C 


T 


C 


0.91 


99.8 


19 


rs7480010 


LOC387761 


11 


42203294 


intergenic 


G 


A 


G 


0.938 


97.4 


20 


rsll 131 32 


EXT2 


11 


44209979 


intron 


G 


C 


C 


0.82 


98.7 


21 


rs 11037909 


EXT2 


11 


44212190 


intron 


C 


T 


T 


0.95 


99.5 


22 


rs3740878 


EXT2 


11 


44214378 


intron 


G 


A 


A 


1 


99.2 


23 


rsl 153188 


DCD 


12 


53385253 


intergenic 


T 


A 


A 


1 


99.5 


24 


rs7961581 


TSPAN8.LGR5 


12 


69949359 


intron 


C 


T 


C 


0.93 


98.2 


25 


rs8050136 


FTO 


16 


52373766 


Intron 


A 


C 


A 


0.24 


99.7 



Chr., chromosome; HW-P, P-value for Hardy-Weinberg equilibrium. 
doi:1 0.1 371 /joumal.pone.0095045.t002 



Weinberg equilibrium (Table 2). The average call rate was 
99.08%. 

Single-locus SNP association analysis 

The results of SNP association analysis with T2DM are listed in 
Table 3. Among all SNPs, six SNPs in five genes were nominally 
associated with T2DM, including rs4402960 in IGFBP2 (OR, 
1.13; 95% CI, 1.033-1.257, P= 0.013), rsl470579 in IGFBP2 
(OR, 1.1; 95%CI, 1.006-1.220, P= 0.049), rsl0946398 in 
CDKAL1 (OR, 1.14; 1; 95% CI, 050-1.249; P= 0.005), 
rsl3266634 in SLC30A8 (OR, 1.22; 95% CI, 1.127-1.334; 
P-CO.OOOl), rsl 11 1875 in HHEX (OR, 1.19; 95% CI, 1.088- 
1.303; P<0.0001), and rsl 153188 in DCD (OR, 1.55; 95% CI, 
1.107—2.18, P=0.02) genes. The directions of association are 
consistent with those reported in European GWAS [3-8] except 
for rsl 153188 in the DCD gene. The reported risk A allele at 
rsl 153188 in European population was associated with reduced 
T2DM risk in our study. The numbers of risk alleles carried by 
T2DM cases were significantly more than those carried by the 
controls (20.94 vs. 20.35 alleles, P<0.0001) (Figure 1). 

Genotype score and diabetes risk 

To examine the cumulative effect of risk variant on diabetes 
risk, we generated weighted genotype score for each risk allele 



from logistic regression model. All participants were divided into 
four equal groups according to genotype scores (<29, 29-31, 32- 
34, and >34). The corresponding OR for each group compared 
with the lowest quartile groups are 1.26 (95% CI, 1.028-1.54, 
P=0.02), 1.31 (95% CI, 1.072-1.61, P= 0.006), and 2.22 (95% 
CI, 1.81-2.73, P<0.001), respectively (Figure 2). We then tested 
whether the incorporation of genotype score improve the 
prediction of T2DM risk. The C statistics of the regression model 
using clinical factors including age, sex, and BMI is 0.627 (95% 
CI, 0.606-0.650). The incorporation of genotype score increased 
the C statistics to 0.657 (95% CI, 0.636-0.678) (P<0.0001) 
(Figure 3). 

Gene-gene interactions and interactions between 
genetic variants and other known T2DM risk factors 

We next explored potential gene-gene interaction and interac- 
tion with other risk factors of T2DM. No significant gene-gene 
interaction was found using pairwise interaction testing. Analyses 
of interaction between genetic variants and other known T2DM 
risk factors showed significant interaction between rsl 3266634 in 
SLC30A8 gene and age on T2DM risk (Pfor interaction<0.0001, 
adjusted P<0.0001). As shown in Figure 4A, the OR associated 
with the risk C allele was attenuated with advanced age, ranging 
from 1.73 in the group of lowest age quartile (age<48) to 0.88 in 
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Table 3. SNP association with type 2 diabetes. 



Power (%) 



# 


SNP 


Gene 


MAF Cases Control 


OR (adjusted 
OR)* 


95% CI (adjusted 95% CI)* 


P /adjusted P)* 


(adjusted 
Power) * 


1 


rs2641 348 


ADAM30 


0.0306 


0.0297 


1.018 (1.029) 


0.757-1.368 (0.762-1.389) 


0.91 (0.85) 


5.19 (5.49) 


2 


rs10923931 


NOTCH2 


0.0302 


0.0288 


1.040 (1.053) 


0.772-1.402 (0.778-1.426) 


0.79 (0.73) 


5.91 (6.59) 


3 


rs7578597 


THADA 


0.0033 


0.0066 


2.00 (1.618) 


0.933-4.29 (0.745-3.514) 


0.075 (0.22) 


80.9 (46.4) 


4 


rs 10490072 


BCL11A 


0.001 


0.0013 


1.346 (1.695) 


0.300-6.024 (0.375-7.656) 


0.69 (0.49) 


7.80 (15.2) 


5 


rs4402960 


IGF2BP2 


0.2551 


0.2267 


1.168 (1.166) 


1.036-1.316 (1.033-1.317) 


O.OII (0.013) 


80.4 (79.5) 


6 


rsl 470579 


IGF2BP2 


0.261 1 


0.2385 


1.128 (1.127) 


1.002-1.269 (1.000-1.271) 


0.046 (0.049) 


60.1 (59.3) 


7 


rsl 70361 01 


SYN2,PPARG 


0.0309 


0.0304 


0.983 (0.976) 


0.757-1.366 (0.724-1.318) 


0.91 (0.87) 


5.17 (5.02) 


8 


rs4607103 


ADAMTS9 


0.3401 


0.3334 


0.970 (0.975) 


0.870-1.081 (0.873-1.089) 


0.59 (0.65) 


9.29 (7.87) 


9 


rs10010131 


WFS1 


0.0599 


0.0613 


1.022 (1.014) 


0.831-1.26 (0.822-1.252) 


0.83 (0.89) 


5.57 (5.23) 


10 


rs6446482 


WFS1 


0.0664 


0.0694 


1.047 (1.038) 


0.859-1.277 (0.849-1.269) 


0.65 (0.71) 


7.91 (6.40) 


11 


rsl 0946398 


CDKAL1 


0.3908 


0.3527 


1.180 (1.167) 


1.061-1.312 (1.048-1.299) 


0.002 (0.005) 


92.2 (88.3) 


12 


rs9472138 


VEGFA 


0.1377 


0.1337 


0.968 (0.968) 


0.834-1.123 (0.832-1.126) 


0.66 (0.67) 


7.59 (7.59) 


13 


rs864745 


JAZF1 


0.2111 


0.2237 


1.074 (1.084) 


0.950-1.212 (0.958-1.227) 


0.25 (0.2) 


24.3 (30.2) 


14 


rsl 3266634 


SLC30AS 


0.4257 


0.4924 


1.310 (1.302) 


1.182-1.451 (1.173-1.446) 


<0.0001 
(< 0.0001) 


99.9 (96.8) 


15 


rsl 1 1 1875 


HHEX 


0.329 


0.2895 


1.199 (1.213) 


1.075-1.337 (1.085-1.355) 


<0.0001 
KO.OOOI) 


94.1 (96.6) 


16 


rs7923837 


HHEX 


0.2054 


0.188 


1.113 (1.126) 


0.980-1.264 (0.990-1.282) 


0.098 (0.071) 


43.96 (50.8) 


17 


rs7903146 


TCF7L2 


0.0234 


0.0287 


0.81 (0.89) 


0.500-1.310 (0.55-1.46) 


0.36 (0.66) 


98.3 (60.5) 


18 


rs290487 


TCF7L2 


0.3965 


0.3810 


1.069 (1.076) 


0.962-1.187 (0.967-1.197) 


0.22 (0.18) 


27.4 (32.6) 


19 


rs7480010 


LOC387761 


0.2196 


0.2163 


1.018 (1.024) 


0.901-1.15 (0.904-1.159) 


0.77 (0.71) 


6.12 (6.98) 


20 


rsl 1 13132 


EXT2 


0.411 


0.4178 


0.972 (0.979) 


0.876-1.078 (0.881-1.089) 


0.59 (0.69) 


8.93 (7.24) 


21 


rsl 1037909 


EXT2 


0.4178 


0.4252 


1.031 (1.022) 


0.874-1.075 (0.920-1.135) 


0.56 (0.68) 


9.67 (7.36) 


22 


rs3 740878 


EXT2 


0.4183 


0.4200 


1.007 (1.001) 


0.895-1.101 (0.900-1.112) 


0.89 (0.99) 


5.24 (5.05) 


23 


rsl 153188 


DCD 


0.0193 


0.0113 


0.583 (0.602) 


0.380-0.892 (0.392-0.925) 


0.013 (0.02) 


79.1 (73.3) 


24 


rs7961581 


TSPAN8XGR5 


0.231 1 


0.2332 


0.987 (0.966) 


0.874-1.114 (0.854-1.093) 


0.83 (0.58) 


5.62 (9.48) 


25 


rs8050136 


FTO 


0.1268 


0.1244 


1.024 (1.018) 


0.877-1.195 (0.869-1.192) 


0.76 (0.82) 


6.27 (5.72) 



MAF, minor allele frequency; OR: odds ratio; 95% CI, 95% confidence interval; * Associated with risk alleles adjusted for age, sex, and BMI. 
doi:1 0.1 371 /journal.pone.0095045.t003 



the group of highest age quartile (age>69). A suggestive 
interaction between BMI and rsl 0946398 in CDKAL1 gene (P 
for interaction = 0.006, adjusted P= 0.084) was also found. The 
increased risk associated with the C allele at rsl 0946398 was 
attenuated in subjects with larger BMI (Figure 4B). However, we 
did not observed significant interactions between genetic score and 
other T2DM risk factors including age, sex, and BMI. 



O Control 
wm Type 2 diabetes 




1112131415161718192021222324252627282930313233 
Number of risk alleles 

Figure 1. Distribution of risk allele frequencies in type 2 
diabetic cases and controls. 

doi:1 0.1 371 /journal.pone.0095045.g001 



Meta-analysis 

Since the statistical power of our current study was not adequate 
for rare alleles or alleles with small effect size, a meta-analysis 

Pfor trend< 0.001 



o 

s? 

m 
2. 
O 

+J 

RS 
i_ 

tn 

T3 
T5 

O 



. I I 



— i — 
<29 



29-31 32-34 >34 genotype scores 



Figure 2. Odds ratio of type 2 diabetes according to weighted 
genetic scores. 

doi:10.1371/joumal.pone.0095045.g002 
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100% - Specificity 0 /. 



Figure 3. Receiver-operating characteristic curve using clinical factors (age, sex, and body mass index) (black line) or clinical factors 
plus genetic scores (gray line) for discrimination of type 2 diabetes. 

doi:1 0.1 371 /journal.pone.0095045.g003 



pooling 20 studies in Han Chinese was conducted to further 
validate the association [13,16-34]. As shown in Table 4 and 
Figure SI, 10 genetic variants in 8 genes are significantly 
associated with T2DM in the meta-analysis, including rs4402960 
and rs 1470579 in IGF2BP2, rs 10946398 in CDKAL1, rs864745 in 
JAZF1, rs!3266634 in SCL30A8, rsl 1 1 1875 and rs7923837 in 
HHEX, rs7903146 in TCF7L2, rsl 113132 in EXT2, and 
rs8050136 in FTO gene. The OR conferred by each risk variant 
did not differ significandy between Han Chinese and European 
populations (Figure S2A). However, substantial heterogeneity of 
risk allele frequencies exists between Han Chinese and European 
populations (Figure S2B). 

Discussion 

In this study, we confirmed the association of 10 genetic risk 
variants identified from European GWAS in Han Chinese. The 
incorporation of genetic information improves the prediction of 
T2DM. The effect sizes conferred by risk variants are similar but 



the allele frequencies differ substantially between Han Chinese 
and European populations. 

Previous GWAS in Han Chinese identified several candidate 
variants associated with T2DM. Tsai et al reported genetic 
variants in PRPRD and SRR genes associated with T2DM in a 
GWAS conducted in Han Chinese in Taiwan [38]. Another 
GWAS in Han Chinese by Shu et al found that genetic variants 
near CDC123/ CAMD1A, SPRY2, and C2CD4B genes are associ- 
ated with T2DM[28]. However, the loci discovered from both 
GWAS did not overlap with each other. Therefore, instead of 
testing these variants, we attempted to validate the association of 
established risk loci in Europeans in our population. Given the 
heterogeneous genetic structure of different ethnic populations, it 
is necessary to validate the relative contribution of T2DM variants 
identified from Caucasian GWAS in Han Chinese. Here, we 
confirmed the association of genetic variants in the IGF2BP2, 
CDKAL1, JAZF1, SCL30A8, HHEX, TCF7L2, EXT2, and FTO 
gene with T2DM in Han Chinese population. Interestingly, the 
effect sizes conferred by these variants were similar between Han 
Chinese and European populations despite marked differences in 
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Table 4. Meta-analysis for SNPs at twenty-four genetic loci in the European and Asian populations. 



Han Chinese Reported associations in Europeans 



Control Power Control 

MAF OR (95% CI) P N (%) MAF OR (95% CI) P N Reference 



rs2641 348 


ADAM30 


0.0277 


1.055 (0.859-1.294) 


0.61 


6,635 


8.6 


0.107 


1.10 (1.06-1.15) 


4.00E-07 


60,048 


8 


rsl 0923931 


NOTCH2 


0.029 


1.027 (0.826-1.279) 


0.81 


5,747 


5.7 


0.106 


1.13 (1.08-1.17) 


4.10E-08 


58,667 


8 


rs7578597 


THADA 


0.0046 


1.300 (0.814-2.078) 


0.27 


6,658 


22.0 


0.098 


1.15 (1.10-1.20) 


1.10E-09 


60,832 


8 


rsl 0490072 


BCL11A 












0.276 


1.05 (1.03-1.08) 


1 -00E-04 


59,682 


8 


rs4402960 


IGF2BP2 


0.262 


1.132 (1.09-1.176) 


5.50E-09 


30,432 


100 


0.29 


1.14 (1.11-1.18) 


8.90E-1 6 


32,554 


4 


rsl 470579 


IGF2BP2 


0.276 


1.109 (1.048-1.174) 


3.49E-05 


13,425 


98.5 


0.3 


1.17 (1.11-1.23) 


1 -30E-09 


13,781 


4 


rsl 70361 01 


SYN2,PPARG 












0.073 


1.15 (1.10-1.21) 


2.00E-07 


59,682 


8 


rs4607103 


ADAMTS9 


0.363 


0.993(0.948-1.041) 


0.78 


15,389 


6.1 


0.239 


1.09 (1.06-1.12) 


1 -20E-08 


62,387 


8 


rs10010131 


WFS1 


0.051 


1.082 (0.952-1.231) 


0.23 


9,370 


25.2 


0.4 


1.11 (1.07-1.16) 


1 .40E-07 


20,922 


6 


rs6446482 


WFS1 


0.051 


1.135 (0.958-1.344) 


0.14 


4,955 


34.1 


0.41 


1.11 (1.06-1.15) 


3.40E-07 


20,922 


6 


rsl 0946398 


CDKAL1 


0.448 


1.206 (1.157-1.257) 


<1.0E-10 


1 9,705 


100 


0.319 


1.12 (1.08-1.16) 


4.10E-11 


32,554 


7 


rs9472138 


VEGFA 


0.125 


1.012 (0.905-1.109) 


0.97 


6,892 


6.91 


0.282 


1.06 (1.04-1.09) 


4.00E-06 


63,537 


8 


rs864745 


JAZF1 


0.272 


1.092 (1.038-1.149) 


6.98E-04 


1 6,996 


97.0 


0.499 


1.10 (1.07-1.13) 


5.00E-14 


59,617 


8 


rsl 3266634 


SLC30AS 


0.421 


1.160 (1.119-1.202) 


<1E-10 


26,074 


100 


0.306 


1.12 (1.07-1.16) 


5.30E-08 


32,554 


7 


rsl 11 1875 


HHEX 


0.308 


1.158 (1.119-1.198) 


<1E-10 


39,986 


100 


0.47 


1.13 (108-1.17) 


5.70E-10 


32,554 


4 


rs7923837 


HHEX 


0.232 


1.176 (1.125-1.230) 


<1E-10 


23,051 


100 


0.377 


1.20 (1.10-1.30) 


7.50E-06 


5,310 


5 


rs7903146 


TCF7L2 


0.0448 


1.454 (1.258-1.680) 


1 .29E-06 


10,391 


100 


0.26 


1.37(1.31-1.43) 


1 .00E-48 


32,554 


4 


rs7480010 


LOC387761 


0.217 


0.979 (0.908-1.055) 


0.58 


9,007 


10.0 


0.301 


1.17 (1.08-1.28) 


1 .20E-05 


5,480 


5 


rsl 1131 32 


EXT2 


0.404 


1.082(1.017-1.115) 


0.013 


9,365 


80.4 


0.267 


1.17 (1.07-1.28) 


2.90E-05 


5,404 


5 


rsl 1037909 


EXT2 


0.401 


1.026 (0.965-1.092) 


0.41 


9,330 


15.1 


0.271 


1.18 (1.08-1.29) 


1 -30E-05 


5,379 


5 


rs3740878 


EXT2 


0.417 


1.057(0.994-1.123) 


0.078 


9,345 


51.3 


0.272 


1.17 (1.07-1.28) 


1 -30E-05 


5,404 


5 


rsl 153188 


DCD 












0.267 


1.08 (1.05-1.11) 


1 .80E-07 


62,301 


8 


rs7961581 


TSPAN8.LGR5 


0.219 


1.044 (0.989-1.103) 


0.12 


15,487 


39.6 


0.269 


1.09 (1.06-1.12) 


1.10E-09 


62,301 


8 


rs8050136 


FTO 


0.134 


1.126 (1.055-1.201) 


3.24E-05 


1 6,690 


98.3 


0.398 


1.17 (1.12-1.22) 


1.30E-12 


32,554 


7 



MAF, minor allele frequency; OR, odds ratio; 95% CI, 95% confidence interval. 
doi:1 0.1 371 /journal.pone.0095045.t004 



allele frequencies, suggesting that the biological actions of these 
variants are the same across different ethnic groups. The effect size 
of several uncommon or rare variants, including rsl0490072 in the 
BCL11A gene, rs7578597 in the THADA gene, and rs7903146 in 
the TCF7L2 gene, are relatively large (OR: 1.34, 1.30, and 1.45, 
respectively). Therefore an aggregate of all SNPs were used for our 
prediction model of T2DM instead of using only common 
variants. 

We found that the addition of genetic information to clinical 
predictors slightly improved the prediction for T2DM in Han 
Chinese. The 3% increment in C-statistics is consistent with those 
observed in Europeans [10,1 1]. Similarly, Xu et al reported that a 
1.6% increase in C-statistics for T2DM prediction in another 
Chinese population using a set of 19 risk variants. Collectively, 
these studies confirmed a small improvement in the prediction of 
T2DM by incorporating genetic information. It should be noted 
that such increment in C-statistics, albeit statistically significant, 
may not be of clinical significance. However, we also found that 
the OR of T2DM in subjects with the highest genetic score 
quartile was 2.22 as compared with those with the lowest genetic 
score quartile. This estimate is in concordance with the 2.60-fold 
increased risk associated with higher genetic scores in a European 
population as reported by Meigs et al. [1 1]. 

An interesting interaction was found between age and an exonic 
variant in the SLC30AS gene. The SLC30A8 gene encodes a zinc 



transporter specifically expressed in pancreatic beta-cells. We 
found that the effect conferred by the risk allele was diminished 
with aging. The underlying mechanism is currently unknown. 
Zinc deficiency has been shown to develop with advanced age 
when the ability to transport zinc is disrupted [39]. Therefore, the 
reduced zinc transporter capability associated with aging may 
mask the genetic effect of &LC30A8 mutation. However, further 
replication is needed to verify this observation. 

Our study has both strengths and limitations which need to be 
addressed. First, this study provides the largest and the most 
updated meta-analysis for T2DM genetic association in Han 
Chinese population. However, this study are still insufficiendy 
powered to validate the association of variants in ADAM30, 
NOTCH2, THADA, ADAMTS9, WFS1, VEGFA, LOC387761, and 
TSPAN8/LGR5 genes with T2DM, probably owing to their low 
allele frequencies and small effects in Han Chinese. Second, 
variants identified by recent GWAS in East Asians were not 
genotyped [40-43]. With the rapidly expanding knowledge for 
T2DM genetics, further incorporation of new genetic variants is 
warranted for in order to enhance prediction. Third, this study 
could not provide accurate estimation of disease incidence because 
of the case-control design. Therefore, the net improvement in re- 
classification could not be estimated. 

In summary, this study affirmed the association of 10 genetic 
loci with T2DM in Han Chinese. Carriers with higher genetic risk 
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scores have a 2.2-fold increase in T2DM risk and the addition of 
genetic information to clinical factors lead to a ~3% increment in 
the discriminative power for prediction of T2DM. These data, 
together with previous studies, support the usefulness of genetic 
testing for T2DM prediction. 

Supporting Information 

Figure SI Forest plots for meta-analyses showing odds 
ratios of type 2 diabetes conferred by risk variants 
identified from European genome-wide association 
studies in Han Chinese. 

(DOC) 

Figure S2 Comparison of odds ratio associated with 
risk alleles (A) or minor allele frequencies (B) between 
Chinese and European populations. 

(TIF) 
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